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Abstract 

Follow-on missions to provide continuity in the observation of the sea surface topography once the successful 
TOPEX/POSEIDON (T/P) oceanographic satellite mission has ended are discussed. Candidates include orbits 
which follow the ground tracks of T/P, GEOSAT or ERS-1. The T/P precision ephemerides, estimated to be near 
3 cm root-mean-square, demonstrate the radial orbit accuracy that can be achieved at 1300 km altitude. However, 
the radial orbit accuracy which can be achieved for a mission at the 800 km altitudes of GEOSAT and ERS-1 has 
not been established, and achieving an accuracy commensurate with T/P will pose a great challenge. This 
investigation focuses on the radial orbit accuracy that can be achieved for a mission in the GEOSAT orbit. 
Emphasis is given to characterizing the effects of force model errors on the estimated radial orbit accuracy, 
particularly those due to gravity and drag. The importance of global, continuous tracking of the satellite for 
reduction in these sources of orbit error is demonstrated with simulated GPS tracking data. A gravity tuning 
experiment is carried out to show how the effects of gravity error may be reduced. Assuming a GPS flight 
receiver with a full-sky tracking capability, the simulation results indicate that a 5 cm radial orbit accuracy for an 
altimeter satellite in GEOSAT orbit should be achievable during low-drag atmospheric conditions and after an 
acceptable timing of the gravity model. 


Introduction 

The very successful TOPEX/POSEIDON (T/P) oceanographic satellite mission has demonstrated the ability to 
monitor the Earth's sea surface topography from space with very high accuracy (Fu et al., 1994). The radial orbit 
accuracy of the T/P precise orbit ephemeris is estimated to be 3-4 cm root-mean-square (rms), which is more than 
a factor of three better than the mission objective (Tapley et al., 1994, Nerem et^l., 1994). With the end of the T/P 
mission in 1997, continuous, high precision observations of the sea surface will be interrupted. Studies of long 
term variations in the sea surface topography, particularly those related to global change, will thus be negatively 
impacted. As a result, the oceanographic community has recommended a follow-on mission to provide 
continuity in the observation of the sea surface topography. 

To link new altimeter observations to existing data sets, it is useful to select an orbit which follows the ground 
track of a previous altimeter mission, such as GEOSAT, ERS-1 or T/P. If continuous, high accuracy monitoring of 
the sea surface is desired, then a logical choice would be the T/P orbit. The unprecedented accuracy of the T/P 
orbit ephemeris is due to the extensive force model development efforts, the robust and precise tracking provided 
by satellite laser ranging (SLR) and the French doppler tracking system DORIS (Nouel et al., 1988), and the low 
atmospheric drag associated with the 1300 km altitude. Combined with the 2 cm precision of the TOPEX 
altimeter, it is an ideal platform from which to study the temporal evolution of the ocean circulation. A follow-on 
mission in the T/P orbit would provide the long term, contiguous observations needed to study such phenomena 
as the secular increase in mean sea level. To date, over fifty ten-day repeat cycles of TOPEX altimeter data have 
been collected, and by the projected mission end the data set will span five years. A follow-on mission in the T/P 
orbit would expand this time frame, making it the preeminent altimeter data set for oceanographic studies. 

Another mission option under consideration is the GEOSAT orbit. With a 108° orbit inclination, a greater 
altimeter coverage of the polar regions is provided. This would allow more complete study of the circumpolar 
currents and the ice fields. Although the GEOSAT altimeter data set is of lower quality than the T/P data set, it 
has been studied extensively and can still provide many opportunities for long term studies. The 800 km altitude 
of the GEOSAT orbit, however, results in increased drag and geopotential perturbations on the satellite, making it 
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. - n an nrbit with the T/p precision. The orbit errors associated with the geopotential can be 
reduced to some extent by tuning the gravhy field with tracking data from the first few months of the mission 
bn S orbi^oi due to atmospheric drag will be more difficult to overcome. Such effects wiU be exacerbated 
biMh'e expected increase in soUr aclivily in the late 1990 s. The ERM orbit would have sunder modeltng 
problems, with possible additional difficulties due to the sun synchronous orbit. 

To Ulus, rate some of the problems already entered with t 

J .. , r p<;iiUs P from the noise level of the doppler tracking systems supporting the satellite (primarily U.S. 
Navv C^NET^mcking plus some French, Canadian and Belgian TRANET tracking). These systems have 

are an order of magnitude large, 

T/P and their measurement noise contributes to the uncertainty in the GEOSAT orbit Similarly, the estimated 
md^^cm rncv of toe ERS-1 orbit is on toe order of 10 cm rms (Korel et al., 1994). Due to the failure of 
on eS only a limited amount of SLR tracking data is available for toe ERS-1 orbit computation and 
efforts are being made to supplement toe SLR tracking with dual satellite cross-over measurements that take 

advantage of toe precise T/P 

pSL°l^XSbuted, S h^h precision observations of toe satellite moHon be available to support future 
missions with low altitude orbits such as GEOSAT and ERS-1. 

The T/P experience has shown that tracking data from a combination of SLR and DORIS or 
Positioning System (GPS) can be used to obtain independent orbits that agree with each other at the 2-3 cm level. 
To illustrate the success of the T/P orbit determination effort, a revised orbit error budget is provided m Table . 
^is^ol^rved'tfia^ A^pre-launch sources of orbit error have been substantially reduced, resulting in an overall 

than three times lower than that of GEOSAT or ERS-1. If a follow-on mission 

occupiestile T/P orbit, it is very likely that this level of orbit accuracy can be maintained using ^y^bmahon 

srsr 

£e E^ol^citainty would be decreased significantly if more SLR tracking were avaikble, but it is the 

cpoT 2 satellite where it is estimated that a radial orbit accuracy of approximately 10 cm rms is attamable 
£32 al 1993) Since DORIS provides only range-rate information, it is difficult to calibrate ^ accuracy 

receiver hardware and software to toe ueartomi wiU make if » even more 
attractive system for the support of future altimeter missions. Projecting the DORIS and the GPS orbit cy 

capability downwards to the GEOSAT altitude is thus an area of interest. 

toe 3 cm rms tewL Clearly su p ( oUo „ton mission to an 800 km alhtude orbit, however, has no, 

repea ‘^ ^ 

accuracy attamable will tend to be domma y on the estimate d radial orbit accuracy, 

extermination strategies are presented fo, mitigating toe effects of drag errors, and a gravity tuning experiment is 
carried out to show how the effects of gravity error may be reduced. 
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Precise Orbit Determination Considerations 


KstoricaU^ uncertamty in the geopotential has been the largest source of orbit error for oceanographic satellite 
nu^ons. However, tremendous improvements made specifically for the T/P mission in modeling the static and 

tTlrlv^ 8 COmP T, S ° f * gravity field have substantially reduced these errors. As a result, orbit errors due 
a 5 e rOU S hl Y equivalent to those due to surface forces at T/P altitude. Such improvement is best 

R^Stefor £e°T7 8 p V™ predicted b y ** 8 ravit Y fie,d covariance (Rosborough, 1986). 

Resute for the T/P pre-launch gravity models, TEG-2B (Tapley et al., 1991) and JGM-1 (Nerem et al„ 1994) and 

Sbk2 P T°t r d f r de S : ^*I' 2 (Nerem et al ' 1994) 3nd JGM ' 3 et a1 ' 1994) are provided in 

Tabk 2_ It is observed that the radial orbit error due to the JGM-2 gravity model, which is the reference model for 

mode^ wh^?H°H b i Pr °f dUCed by ^ ASA/ 15 ab ° Ut 2 Cm ' 7116 JGM ' 3 model ' which 15 a refinement of the JGM-1 
l f l ^formation from tracking of T/P (SLR, DORIS, and GPS), Stella, Lageos-2, and SPOT-2 

GEQSATorh > Cm T ° rblt JGM ' 3 m ° del ako y ieIds improvement for the ERS-1 and 

GEOS AT orbits, with predicted rms radial orbit errors of about 4-5 cm. Such will be the level of pre-launch radial 

S ° kl y to , the geopotential for a follow-on mission that occupies either of these orbits. Previous 

tL^Lr? e S ° luh T S ha ^ e dem ° n strated that post-launch tuning of the gravity field can reduce this error, but 

Sfecte ' i? JT 'IT f P r °" l a ^ Uality and ^ abiHty t0 Separate S ravit y effe cts from surface force 
effects. It should also be noted that additional radial orbit error can be expected from the time-varvine 

component of the gravity field. Bettadpur and Eanes (1994) have shown that at the 1300 km T/P altitude, therms 

“ lal orbit err or due to ocean tides is about 1 cm, and increases to about 3 cm at 800 km altitude. These errors 

ay be reduced through improvements of the ocean tide model with altimeter data being obtained from T/P. 

ha V a alS ° bCen la u rge sources of orbit error for oceanographic satellite missions, especially 
the 1 33flT df |H 8 H 7 ma) ° r T tribUt ° r t0 ^ T/P ° rbit P recision » low atmospheric density associated with 
mfJnn 0 ™ 3 ^ f ^ ^ *7 Telativel y low soIar and geomagnetic activity occurring throughout the 
mission. These elements have combined to provide T/P with a low drag environment favorable for computing 

r xam r'% at ,' he T/F alHtUde ' ,he abnos P h " e tenuous. and the accelerations duf 
A ^ d thermaI mdiation pressure are one to two orders of magnitude greater (Ries et al 19921 

At the 800 km altitude of GEOSAT and ERS-1, where the atmosphere is much nTe dens^ drag aSeLmtio^ may 

bX^T 6 40 ra i a 5 0n P T, SUre ' T ° Ulustrate 411656 Concepts ' t y pical accelerations due to surface forces at 

each month kiT992TnH de<1 h ^ accel crahons were computed from daily averages for the first day of 

each month m 1992 and are based on the precise force models used in the UT/CSR orbit determination software 

terresto.7md*f * ** 13 °° ^ alhbld \ the acceleration due to drag is much smaller than that due to solar and 
terrestrial radiation pressure, and that the variation in the drag acceleration from minimum to maximum is 

usually no more than an order of magnitude. However, at the 800 km altitude, the acceleration due to drag can 
f^ 666d ^ at d “ e t0 Solar ^rrestnal radiation pressure and exhibits very large fluctuations with the maximum 
being two orders of magnitude greater than the minimum. Despite the fact that the acceleration due to drag is 

deStTfl h ^ 7 f UC t0 radiation pressure ' if 15 dra 8 *at has the most dramatic effect on the orbit. Rapid 

^ ° CCUr SeV6ral 3 day makC m0d6lin S effects of Aspheric drag very 


Simulation Description 

Hie simulation performed in this analysis is intended to provide a preliminary assessment of the radial orbit 
accuracy that can be achieved for a satellite in GEOSAT orbit, with emphasis on the impact of gravity and drag 
d7nend nd t h ° W ^ s y ste m technology can be used to overcome them. The fidelity of the simulation is 

thfi rtf UP ° n r6a 7 hC S ° UrC6S ° f d y namic and measurement model errors, and effort must be made to ensure 

rn^kh llTf? 3re 116141161 T rly optimistic nor Pessimistic. Thus, it is necessary that the orbit error 
™ d6ls , b cahbrated against those observed in real data analysis. The dynamic and measurement model errors 
used m this simulation are generated through differences between the models used in the data generation and the 
ata processing. Random stochastic effects are added in the data generation phase. The residual rms radial orbit 
ference that remains after the data is processed is the radial orbit error. This simulation makes extensive use of 
the error models developed for GPS applications by Rim et al. (1993). Improvements to some of these error 
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models were made with updated Mormation from 

the GEOSAT orbit were developed for this analysis with knowledge gamed trom me proc g 
data, especially T/P and SPOT-2. 

Calibration of the dynamic error model is carried out by 

reference force model and then processing e resu g ft position or "idealized" tracking, and thus 

equivalent to having continuous, perfect observaho^ of t^ spacea^ p^ ^ si ^ ulated b gener ating a satellite 

considers only force model erTOTS ' ° f * e n processing the resulting satellite positions with the JGM-2 

ephemeris using the JGM-1 gravity model and tten processugme res g ^ a ^ articular source of error, 

model. Adjusting only the satellite ep s a e pr IGM-1 and JGM-2 Adjusting additional parameters 

which in «. ~ in « or accommoda^g 

SS™rr™"« considered. Most of to measurement mode, errors were 
generated based on uncertainties obtamed from the current literature. 

The GPS measurement assumed for the simulation was 

double-differenced mode. Effects of selective aval a tton / /or ^ (light receiver and the ground network 

%%%£££ ’to Table' T pTrSom^was assumed that a GPS flight receiver capable of tracking all GPS 
satellites in view will be available to support the orbit determination. 

An assumption on to satellite design is also neces»ry 

to-mass ratio and solar array size are necessary to properly scale the T/p> sPOT-2 and 

descriptions of some relevant satelhtes are provi e m spites are configured in a "box-and-wing" design, 
designs are similar in maximum area-to-mass ratio Tb«e is* teflttes are co, itigu o( a box with „ attached 

where the spacecraft is heated as a com ma ton altitude option because its satellite bus is very similar 

sour array. The SPOT-2 model was chosen for the ! 800 to < altttade °P«° sma ller solar , is als0 

as ERS-1. 


Error Model Description 

The important sources of dynamic model error for r parti^lariylmportot. Gravitational errors 

forces, and at GEOSAT altl ^ e ' * e ^arth^gravifational coefficient, GM, errors in the geopotential and ocean 
included a bias m the knowledge ° & , , Farth tide parameter. The errors for the geopotential are 

tide model coefficient and a bias m tiie dymanuc ^sohd Ea * 0C ean tide m J el are bas ed on 

based on JGM-2 and TEG-2B gravity model differences, ana tne errors 

differences between selected ocean tide models. 

Non-gravitational errors for both the altimeter satellite 

“enormode, employed /tills simulation is summarized UTahle 

The important sources of observational error ^e ins.^ mSur^en. 

tropospheric refraction and ionospheric disp _ narameters nutation and precession errors, and station 

OTwdinate errore.* I ^^basic < ^t^ type^Tsumed was JgPS dual-frequency carrier phase measurement in the 
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double-differenced mode. The effects of clock error including selective availabilitv mulH nat+i 
migrahon and ionospheric dispersion were neglected. Table 7 summarizes the measurement model erreS CG ° 


Error Model Calibration 

As previously discussed, the dynamic error model consists of force model differences between the perturbed orbit 
and the reference orbit, with random stochastic effects added. The level of orbit error generated in this manner is 
^ Tiff 31 ? 8 the s P ac , ecraft e q u a«ons of motion with the reference force model and then processing 

S? ♦ ? W u u / ertUrbed m ° deL °y namic model Parameters are estimated over various arc 
d sub-arc lengths to absorb the force model errors, and the residual radial orbit differences that remain are 
compared to the expected level of orbit errors derived from real data experience. 

^SinZ H eSU n S the C l ynamiC erf0r model 15 8 iven Table 8, where the dominant error sources are 
examined individually. For each case considered, a ten day ephemeris was generated with five minute soacine 

:“ POd *K Were utiHzed to “SW* effects of low and high atmospheric activity on the drag SS 
model. Hie ephemens was processed such that the satellite state was adjusted once over the entire ten day arc 

^T empmCaI ir? r t f.r SVerSe (T) 311(1 n0rmal (N) Parameters were adjusted for each day. In the case of the 

m °fr r bra T' ! 6 P arameters were adjusted every six hours, while for T/P, along- 

tr^ck accelerations (Cy) were adjusted every six hours. , ® 

* e e " or due to J GM ' 2 TEG-2B differences is examined. The model differences for these 

£ <£re A TaZl A ^ *° * T/P alB,ude - 1>P“ '» -ry reaWc radiXSS 

nn anri a *J tad Table 2 ' ^ ^ radia l orbit error predicted from the JGM-2 covariance is 7 4 

ueTlcM^d IGmTh^ B ined Wilh * he modd differenres yid* the result In Case 2, the gravit^ 

"Z°Z T/ J n r ^ d 11 dd ' ere " CeS 15 c * ami ™ d model deferences result in 2.1 cm rms radial orbit 
/ , which is compared to value of 2.2 cm from the JGM-2 covariance. The 3 6 cm rms radial orbit 

covarianr nerat H d ^ GEOSAT with *** model deferences is about half of that predicted by the JGM-2 
covariance and is more representative of what tuning the gravity model might provid e P Accordingly the errors 

fio? J d r y * loremendoned model differences will be termed the "pre-launch" and "tuned" gravity model 

Stade 6rr0r ^ r6SUltS ^ 3bOUt 05 0111 ° f radial orbit error re gardless of die ort>it 

^1% ^, C3 4 ' 006311 5ohd Earth fa de errors result in about 3 cm of radial orbit error at GEOSAT altitude 

In Case 5, the drag error model was examined. It is observed that drag is negligible at T/P altitude regardless of 
die epodi. However, at GEOSAT altitude, the rms radial orbit error LreaJdby more than a 

as lar^/s rt! mble adnos P h ® re of Ma rch, 1991 epoch. Moreover, this level of radial orbit error is about twice 
l a H" ^ e , S ray dy error from the "tuned" model and is commensurate with the gravity error from the "pre- 

obS^vedTaflevel ^ ^ 5 e , SOl £'l terrestriaI ' and thermal radiation pressure error model is examined. It is 
observed that level of rms radial orbit error, which was generated from a fairly pessimistic error model is small in 

comparison to both the gravity and the drag error. This again demonstrates the ability of Ae emS 
accelerations to absorb mismodeled or unmodeled accelerations, especially those that have distinct 1-cpr 
signatures such as radiation pressure errors. ^ 

CffeCt ° f 411686 in ^^, dUa ^^ amiC err ° r SOUrCes ° n * e 11115 radial orbit accuracy is shown in Table 
* { m A al dynamiC error model for GEOSAT orbit was selected to emulate "pre-launch" dynamic errors and 

lh e J^H dyna T^ 6 7° r m ° deI f ° r li P ° rbit WaS selected to emulate the current dynamic errSs that have been 
l S 17! d m " ea data P r °° essm 8' Tb e effect arc length was investigated to illustrate the impact of dynamic 
model error build-up on the rms radial orbit accuracy. Long arcs, on the order of a few to several days iSwide 
dynamical constraints through the satellite equations of motion, but suffer from the build-up of non-con^rvative 
surface force model errors. Short arcs, on the order of a day or less, tend to attenuate long P period and re^onanl 
pavity errors and reduce the build-up of non-conservative surface force model errors. For each arc length used 

hours X d ; tl0nS T Wer6 ad i USt6d at . arc e P° ch C D 's for GEOSAT and Ct’s for T/P were adjusted fvery six 
hours. The 1-cpr T and N accelerations were adjusted daily, except in Case 4, where they were adjusted every six 

£n U d S a t f aS6S J° US 3 'u 4 observed that for both T/P ^d GEOSAT shortening the arc length froi^one 
ten-day arc to ten one-day arcs has little effect on the rms radial orbit error, although a slight benefitwas gained 
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for GEOSAT in ft. March. 1991 thne frame wh 

arcs essentially contain all the dynamic orb acce i era tions This implies that long arcs, such as 

of the 1-cpr T and N accelerations in absorbing slowly changmg a manner that does not 

the ten-day arc length used to generate the precise T/Po^™ “ 1 ’"^4 tesubTrc length of the 1-cpr T 
result in significantly increased orbit error !^I% to J^eSonb3S, Approach the "reduced dynamic" 
and N parameters is reduced to six hours, P (Bertieer et al, 1994). It is observed that the 

filtering technique used success u ^ «‘ hG reduction ^ the rms radial orbit error generated for both 

parameterization in Case 3 results m a sign March 1991 time frame. However, only a 

GEOSAT and T/P in the February, 1993 time frame, and for T/ GEQSA t ^ dicating that additional C D 's are 

negligible reduction in the March, 1991tlI ^ e r ^ ™ observed that the time frame of the simulation has little 

- a, ***. . 

one day or longer are used. 


GPS Error Model Calibration 

Calibration of the d,tt S^d^ We^ay pedS 

one utilizing the classical approach and the oth g ^ ^ firsl casc W tial conditions, a Y-bias 

the twelve days of data were processe V , , were adjusted at arc epoch- The real and 

parameter and a scale factor for the box-wing radiation .P^^' ”^“’ r 0 72.2 to Jem. The smallest 
simulated data fits showed excellent agreemen . wi 0 01 cm, and the largest was 0.1 

difference in residual rms observed between Ute real and ^sunulated febfe w£ aoi c ^ „„ 

cm. In the second case, initial conditions, a C T and a 1-cpr T ' and IN M , c [ 6 cm pne smallest 

simulated data fits again showed exce ent a greemen / 0 12 cm In either case, the simulated data did not 

difference observed was 0.02 cm and die largest mSeTaJnot too optimistic. In Tables 10 and 11, 

St consistently below that of the tea data, 1 md'oabng th >bove ustag Jv-bias and radiaHon pressure 

eBecTfortiie^re^ndfonam of d^^uU^nwnuUt^ quite^weU dmtvvt^a^ observed in real data processing. 


simulation Results 

Eighteen days of GPS double-differenced phase data were simu^or 

orbit beginning February 2, 1993 with t e compre ® conditions conducive to the separation of drag effects 

described. The epoch was selected to provi ea P To test the effectiveness of tuning the gravity 

from gravity effects, so the results will be days of the data were obtained with 

field for a satellite in this orbit, information d f h , unt .J field (or comparison purposes, 

four three-day arcs. This prov,ded s« dg ^dam ““ e ^ s ! multaneously alo „ g with die position of 
Orbits for the altimeter satellite and a with the smallest random position errors were held 

twelve of the fifteen GPS ground stations. The six-hour^ and a daily 1-cpr T 

fixed as fiducial sites. The parameterization used for tt* Recall that JGM-2 is the reference 
and N, while that for the GPS sate Utes was a model difference. The recovery 

gravity field and TEG-2B is the field used to g en «ate gravity^ d ^ dif f erence variance (DDV) in Fig. 1. The 
of the JGM-2 field in this experiment is shown yp ^ ^ ^8 tQ JGM . 2/ and the dashed lines with open 
heavy line with solid dots represents the gravity held were perfectly recovered, the DDV would 

be°tre P STd ^ It S that TEG-2B and JGM-2 agree well for the dominant zonal terms such as J* 
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to fte' ZZgpZZ l7dZ™»X%iT2?’Z‘ OUT ' JTt reC0Very °‘ the ,GM - 2 mod '‘ * 
resonant terms. Above degree 24 the DDV of Hip f a r u mucb of the improvement attributed to the 

suiting ft* some of ( °S ^ U " tii ab »“> <«•« % 

eighteen days of data are included in the gravity mod^l thp 7 Contammated the gravity model. When all 
that tuning the gravity field with the satellite tracking Hal* ■ u ^ Were es f :ntia lly the same. This demonstrates 
geopotential, particularly ^ ** radiaI ° rbit error due *> *e 

tracking data were used to tune the graviW fieW fn 1 ^ ^ ° rder terms ‘ ° nly twelve da y* of 

with the additional information provided by a few mtSTh ' n fur * he '_ improvement could be realized 
redundancy required to average ou, the random components of the da,a 

^differenced phase data with the pre-hmed and 
the daily Y-bias and three-dav radiation ^1 p fJt da ys of data were processed m six three-day arcs using 
and daily 1-cpr T and N paraLterizahonfoJ iS P a ; am ««^>hon for the GPS satellites, and the six-hour Cn 

nns the Z 5 7 cmt^ T T r "*“ d ^ * 

rms with the post-tuned field. It is interestine to note that aim * S ,' d * subse< ? uentI y reduced to 3-5 cm 
tune the gravity field, improvement in those residuals was stil/nhi - h ^ tbree ~ days arcs were n »t used to 
orbits to the reference orbit shows that radial orbit errors of about 8-1 1 cm ” erm ° re ' com P arison of the best fit 
gravity field and all other dynamic and measurement errors wMetbou^ ZTJZ? ^ ** P 1 ^^ 

post-tuned field. The reduction in rms radial orbit error due to gm T ***** With 

it is considered that only twelve da vs of a sin vie 17 A™ rpncfx ^ 1S ^ S1 gnificant, especially when 

improvement is obtained ^ ^ ^ ***« 

13. It is observed that for both the pre-tuned and w 3!L antnbution of drag error, as shown in Table 
when one-day arcs were used The reduction in i u *. ' resi duals are reduced roughly by half 

lha, short ani on the ord^of o^ day wT.'f the 5ama ,OT b « h “4 indicating 

GEOSAT altitude. V ‘° con,end w,th the ^el of drag error associated with the 


Conclusions 

GE^^T^rbit^h^^acl^d h^GP^T^^iunwic^^s^uIatiOTi^echni^ ^ ^^a f ° r “ aItimeter SatelUte * 

of orbit error with dynamic and measurement error models ralih a t0 generate reaIisti c sources 

The GPS tracking scenario assumed a constellation of 21 Block II satellites witTfl^hi 1631 < ? t ® P rocessin S- 
can track all GPS satellites in view. A gravity tunine exnerimpru ^ ' T* fl ! ght a ° d ground receivers that 
radial orbit error due to the geopotential can be redureH a* Tf “ nducted to demonstrate how the rms 
examine strategies effective in redu^ 3nd ° rblt determ ^ a tion tests were conducted to 

expected radial orbit error due solely tcfthe geopotential canbf reH^ ^ S , tUdy ' 14 Was found that the 

field. Only twelve of the eighteen dLs of simnK l ! edUCed Wlth P ost - la ^nch tuning of the gravity 
reduction L the radial orbft e^or wa T* *° tUn6 ** yet a signdic ^ 

difficult to mitigate. The simulation ^ FV™ * *° dr3g ' hOVVeVer ' wiI1 be more 

field, about 5 cm of rms radial orbit error can be exnertpH h ° W F j g environment with a tuned gravity 
robust tracking is provided bv a GPS flivhi ^ ui W one ' day arcs are use d- This also assumes that 

error model cSibration showed Zt t llTeZ Zl ° f traC, T 8 aU GPS SateUites ® view. The dynamic 

variations exist. Adjusting the geopotential" mode? SubstantiaUy when hi g h atmospheric density 
topic for further study. P ° del m 11115 hl S h ' dra g environment will be more difficult and is a 
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Table 1. TOPEX/Poseidon Orbit Error Budget 



Table 2. RMS Radial Orbit Error Predicted by Gravity Covariance 


TEG2B 

JGM1 

JGM2 

JGM3 

TOPEX 

(cm) 

10.1 

3.4 1 

2.2 
1.0 

ERS-1 

(cm) 

19.2 

8.3 

7.9 

4.0 

GEOSAT 

(cm) 

15.3 
8.2 
7.4 

5.3 


Table 3. Typical Accelerations Due to Surface Forces 


Satellite 

Altitude 

(km) 

Average 
Solar 
Radiation 
(10” 9 m/sec 2 ) 

Average 
Earth 
Radiation' 
(10“ 9 m/sec 2 ) 

Min / Max 
Atmospheric 
Drag 

(10“ 9 m/sec 2 ) 

T/P 

1300 

! 60 

3 

0.2 / 1.7 

GEOSAT /ERS-1 

780 

70 

9 

3/136 


Table 4. GPS Tracking Assumptions 


• GPS constellation consists of 21 Block II satellites 

Flight and ground receivers track all GPS satellites in view 

15 elevation cut-off for ground stations; 5° elevation cut-off for altimeter 

• 15 GPS ground stations 

• 2 minute observation sampling time 
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Table 5. Altimeter Satellite Orbit and Spacecraft Geometry 


I 

Projecte 

d Areas 


Satellite 

T/P 

cpoT-2 

Alt. 

(km) 

1330 

820 

Inc 

(deg) 

66 

98.7 

■1 

mi\VM 

Max 

Area/mass 

(m 2 /kg) 

0.012 

0.013 

Roll 

(m 2 ) 

4.7 

6.5 

Pitch 

(m 2 ) 

8.2 

3.5 

Yaw 

(m 2 ) 

8.3 

9.0 

Solar 

Panel 

(m 2 ) 

25.5 

18.5 

AO A 

jl vy i a. 

ERS-l 

GEOSAT 

770 

780 

98.5 

108 

2400 

662 

0.014 

0.008 

5.0 

5.0 

5 

15.0 

.0 

28.0 


Table 6. Dynamic Error Models used in Simulation 


• Altimeter satellite and GPS gravitational errors : 

- GM error of 0.0008 km 3 / sec 2 

- Geopotential errors from JGM-2 vs. TEG-2B 

- Ocean tide errors from CSR tide model differences 

- Solid earth tide errors from 3% error in k2 

• GPS non-gravitational errors : 

- Radiation pressure 

> box-wing reflectivity randomly perturbed by 5% 

> thermal acceleration with amplitude of 1 nm/sec , on vs. off 

> 5 km random error in Earth shadow radius 

> integration step size, 500 sec vs. 600 sec 

- Y-bias acceleration stochastically perturbed by 5%, 6-hour correlation time 

- Constant 1° solar panel misalignment, constant 1° solar panel pitch angle error 

• Altimeter satellite non-gravitational errors : 

- Radiation pressure 

> 50% error in solar reflectivity. Earth albedo and emissivity 

> Earth albedo and emissivity also randomly perturbed by 3% 

> 50% error in thermal model acceleration and time decay 

> 5 km random error in Earth shadow radius 

> integration step size, 50 sec vs. 60 sec 

- Drag error from constant flux, random 3-hour K p vs. standard flux and K p 

- Constant 1° solar panel misalignment, constant 1° solar panel pitch angle error 


220 


















Table 7. GPS Measurement Error Model used in Simulation 


• Observational Errors : 

- 5 mm random error for each phase range 

- time tag errors from two-point Allan variance model with 
CTf = 5xl(r 13 , T\ = lxlO 3 , x 2 = lxlO 5 for TPFO receiver 

af = lxlO' 13 , ti = lxlO 5 , t 2 = lxlO 7 for GPS ground receivers 

- troposphere error from modified Hopfield model vs. Chao's model, with 0.8% 
stochastic tropospheric biases added 

- Phase center errors via satellite attitude errors of 1° in roll, pitch, yaw 

• Random precession and nutation errors 

- 0.1 mas/ yr random error in precession 

- 1.0 mas noise in 1-day values of nutation 

- 0.1 mas random errors in long period of nutation 

- 0.05 mas random errors in short period components of nutation 

• 1 mas random noise in 5-day values of Earth orientation (xp, yp, UTl) 

• 3 cm random errors in station coordinates 

• 3% random errors in individual station tide corrections 

• 8 mm/ yr random error in all tectonic plate velocities 


Table 8. Altimeter Satellite Dynamic Error Model Calibration Using 10-Day Arcs 




RMS Radial Orbit Differences 
(cm) 



Epoch 2/1/93 

.. 

Epoch 3/12/91 

Case 

Error Source 

T/P 

Geo 

T/P 

Geo 

1 

"pre-launch" gravity model error 

n/a 

7.4 

n/a 

7.4 

2 

"tuned" gravity model error 

2.1 

3.6 

2.1 

3.7 

3 

GM error 

0.5 

0.5 

0.5 

0.5 

4 

ocean and solid Earth tide error 

1.0 

3.4 | 

1.2 

3.3 

5 

atmospheric drag error 

0.0 

1.8 

0.1 

7.1 

6 

radiation pressure error 

0.2 

0.7 

0.1 

0.5 


Table 9. Final Altimeter Satellite Dynamic Error Model Calibration 



RMS Radial Orbit Differences 
(cm) 



** mm ED cBSl 


Case 

T/P 

Geo 

mam 

Geo 

1 

2.4 

8.9 

2.4 

10.6 

2 

2.4 

8.6 

2.4 

10.1 

3 

2.3 

8.1 

2.3 

9.5 1 

4 

1.6 

6.7 

1.5 

9.4 
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Table 10. Error Model Calibration for GPS Satellites 
3-day vs. 1-day Orbit Solution Differences using Real Data 


GPS 

Satellite 

Radial 

RMS 

(m) 

Transverse 

RMS 

(m) 

Normal 

RMS 

(m) 

RSS 

(m) 

1 

0.2634 

0.7123 

0.3910 

0.8541 

2 

" (14254 

1.1245 

■IK 


3 



0.3004 

0.6563 

11 

0.1980 

0.6501 

0.3796 

0.7784 

12 

0.2324 

0.6680 

0.3091 1 


13 


0.5326 

0.2975 

0.64ZZ 

14 

0.2656 

'fflwrwyi w if 


0.6705 

ft 

“ <53764 | 

—ii 

0.4113 


16 

0.2787 

0.8108 


0.922U 

17 

0.3383 

0.8126 

0.3526 

0.9482 

18 

03738 


0.2630 

0.5096 

19 

0.1763 

0.4771 

0.2488 




1.2239 

0.4165 


21 


0.7410 

0.2823 


23 

0.1971 

“ 53730 



24 


0.7115 

0.3656 


25 ' 

53668 " 

0.8509 

0.3786 

0.9688 

26 

0.2118 

0.6973 

0.3655 

0.8153 

27 

0.2113 

0.6229 

0.3151 


28 

0.2606 




29 

0.2735 

I 0.5825 


1 0.7342 | 


Table 11. Error Model Calibration for GPS Satellites 
3-day vs. 1-day Orbit Solution Differences using Simulated Data 


GPS 

Satellite 

Radial 

RMS 

(m) 

Transverse 

RMS 

(m) 

Normal 

RMS 

(m) 

RSS 

(m) 

1 

0.2021 

0.6088 

0.3536 

0.7325 

2 

0.4454 

1.2733 

0.5832 


3 

“ 0.2526 

0.5433 

0.3302 

kimik 

11 

0.1878 

0.6441 

0.3184 

1 

12 

0.1812 

0.5483 

0.2990 


13 

0.1578 1 

0.4224 

0.2212 

u. 30 zz 

14 

0.2369 

0.4815 


0.5880 

15 

0.2699 

0.7739 



16 ”1 

0.2676 

0.7449 

0.2984 

0.8459 

17 

0.3427 

0.8260 

0.3077 

^ i| 

18 

0.1478 

0.3470 


U.46U1 

19 

0.1744 

0.5760 



20 


1.3548 

0.4643 


21 

0.2638 

0.6340 

0.3163 


23 

0.1737 

0.5282 

0.2835 


24 

0.2604 

0.6162 



25 


0.8529 

0.4322 


26 


0.6022 

0.3692 


27 

0.2017 

0.6303 

0.3061 


28 

0.2911 

0.8480 

0.3687 


29 

0.2651 

0.5459 

0.3443 

i u.by/o | 
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Table 12. Results of GPS Tuning Experiment using 3-day Arcs (solar activity of 2/1/93) 




Pre-Tune 

L Post-Tune 1 

Arc 

1 

Residual 

RMS 

(cm) 

5.7 

RMS 

Orbit Differences 
(cm) 

Residual 

RMS 

(cm) 

RMS 

Orbit Differences 
(cm) 

R 

8.7 

T 

34 i 

N 

17 4 

R 

T 

N 

2 

5.5 

9.8 

29.6 

1/.** 

22.9 

3.6 

4.0 

5.9 

6.2 

17.3 

18.4 


3 

6.2 

10.8 

37.8 

20.6 

4.0 

6.7 

17.1 


4 

7.0 

10.0 

42.8 

21.4 

4.6 

6.5 

19.8 


5 

6.2 

7.5 

40.0 

17.6 

4.5 

5.9 

22.3 


6 

5.5 

7.9 

27.5 

15.5 

4.8 

6.8 

27.4 

R 6 




0.0 


Table 13. Results of GPS Tuning Experiment using 1-day Arcs (solar activity of 2/1/93) 



Pre-Tune 

Post-Tune 1 

Residual 

RMS 

(cm) 

RMS 

Orbit Differen 
(cm) 

ces 

Residual 

RMS 

(cm) 

RMS 

Orbit Differences 
(cm) 

R T 

r N 

R T N 

3.3 

8.2 25.5 

18.5 

1.8 

4-5 12.6 8.6 
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